Conservation of dielectric constant upon amorphization in perovskite oxides 



OO 

o 
o 

(N 



Pietro Delugas,^ Vincenzo Fiorentini,^'^ Alessio Filippetti,^ and Geoffrey Pourtois"^ 
1 ) NXP Semiconductors, Kapeldreef 75, B-3001 Leuven, Belgium 
2) SLACS-CNR-INFM (Sardinian LAboratory for Computational materials Science) and Dipartimento di Fisica, 
Umversitd di Cagliari, Cittadella Universitaria, 1-09042 Monserrato (CA), Italy 
3) IMEC, Kapeldreef 15, B-3001 Leuven, Belgium 
(Dated:) 

We report calculations indicating tfiat amorphous RAO 3 oxides, with R and A trivalent cations, 
have approximately the same static dielectric constant as their perovskite crystal phase. The effect 
is due to the disorder-activated polar response of non-polar crystal modes at low frequency, which 
compensates a moderate but appreciable reduction of the ionic dynamical charges. The dielectric 
response was studied via density-functional perturbation theory. Amorphous samples were generated 
by molecular dynamics melt- and- quench simulations. 
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I. INTRODUCTION 

In the context of integrated microelectronics there is 
a clear dichotomy, rooted in the fundamentals of dielec- 
tric screening, between crystalline and amorphous dielec- 
tric layers. The latter are preferable in terms of electri- 
cal leakage behavior, but their dielectric response may 
be unsatisfactory for applications requiring a large static 
dielectric constant Ks- A case in point is that of the so- 
called high-K materials, gradually replacing silica^ as gate 
dielectric in aggressively downscaled Si-based integrated 
microelectronics. It is unclear in general if the high n of 
a crystal compound will be conserved in its amorphous 
phase. This is due to the nature of the dielectric screen- 
ing in high-K materials, where the total static dielectric 
constant ^s-— ^oo~t~^ion 

is dominatedSii^i^ by the lattice- 
vibrational component Kion, being typically a factor 2 to 
10 larger than the electronic Kqo- Large values of Kjon 
are due to large effective charges and soft infrared (IR) 
vibrational modes. This is borne out by its definition: 
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with z the mode charge vector, VL the system volume, 
Z* ^ p the effective or dynamical or Born charge tensor 
and Mi the mass of atom i, ^^nd uj the eigenvector 
and eigenfrequency of mode A at zero wave- vector, and 
a, /3 cartesian indexes. 

Dynamical charges are termed " anomalous" if they ex- 
ceed the nominal ionicity. Such "anomaly" (sort of a 
misnomer, being ubiquitous and dramatic in many ox- 
ides) results from a subtle interplay of covalent and ionic 
bonding, as discussed in detail in Ref.d. Lattice disor- 
der will generally reduce the efficiency of orbital overlap 
and, hence, the dynamical charge flow measured by the 
dynamical charges. Thus the Z*'s are expected to de- 
crease in the disordered amorphous environment. On the 
other hand, softer IR-active modes counterbalancing this 
effect (through the l/w^ factor in Eq. 1) may appear in 
the amorphous phase (they indeed do, both in the present 



case as well as in zirconia, as reported in Ref.p). The ma- 
jor reason for this is that IR-active and inactive modes 
are not discriminated by symmetry constraints as in the 
crystal phase, so that non-IR crystal modes, or remnants 
thereof, can acquire IR activity. For sufficiently low fre- 
quency, thanks to the l/cj^ factor in Eq. 1, even weakly 
IR-active modes may produce a large dielectric intensity 
(counterbalancing the general decrease in Z*) and cause 
a conservation or a boost of the dielectric constant. 

From the discussion so far, it is apparent that a quan- 
titative assessment of the dielectric constant ingredients 
is central to evaluating a technologically relevant dielec- 
tric. Indeed, if the amorphous phase of a given material 
has a similar ionic screening as its crystal phase, there 
is hardly any point in pursuing the fabrication of single- 
crystal layer of that material. 

Here we study from first principles the static dielectric 
response in amorphous high-K RAO3 oxides with A, R 
trivalent cations. We show that indeed the dielectric con- 
stant in the amorphous phase is similar to that of their 
distorted-perovskite crystal phase, due to a moderate re- 
duction in polarizability, and to the disorder-induced IR 
activation of non-polar low-energy modes of the crystal. 

Our specific case studies are LaAlOs and DyScOs. 
We described previously their crystal properties.^'** The 
nomenclature "R" for the larger, nominally 12-fold coor- 
dinated cation and "A" for the octahedrally-coordinated 
one is used throughout the paper. These two com- 
pounds represent quite thoroughly the relevant trivalent- 
cation perovskites RAO3 for high-K applications. In- 
deed, the choice of A can be restricted to either Al 
or Sc because rare-earths, Lu, La, and Y function as 
the R cation, mixed-valent Ti and Mn cause small gaps 
and correlation)^ and Ga is analogous to Al in this 
stoichiometry.^'^ Both materials are actively studied by 
the high-K community, and both have been reported^ii^ 
to have high dielectric constants (^^25) in the crystalline 
as well as the amorphous state. In DyScOa the effect is 
so far uncontroversial. A report^-'^ also exists of a sizable 
(40%) deterioration of the dielectric constant in LaAlOa. 
While not aiming at a definite resolution of this issue. 
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here we find that a conservation of dielectric screening in 
these materials is plausible from the theoretical point of 
view. 

II. TECHNICAL MATTERS 

Amorphous LaAlOa and DyScOa samples were gen- 
erated by a melt-and-quench molecular dynamics proce- 
dure, using the empirical-potential GULP code.^^ A liq- 
uid was equilibrated at temperature T=5000 K for 30 ps 
to achieve fully randomized cation positions. It was then 
quenched by Nose constant-T dynamics at 50 K/ps. 2-ps 
contant-T intervals separated by abrupt 100 K down- 
steps were used down to 2000 K. After that, since the 
mobility was then very low, larger down-steps of 500 K 
and simulation periods up to 5 ps were used down to 500 
K, whereby structural properties (expected to be well 
converged'''"' at this cooling rate) were sampled in a 5 ps 
run. 

For LaAlOa, we also ran ab initio density- functional 
constant-T dynamics with the VASP code^"* to quench 
the equilibrated liquid from 5000 to 500 K with alternat- 
ing 0.5-ps constant-T intervals and 0.5 ps, 500 K ramp- 
down segments. The differences in structural properties 
in the two approaches are moderate and affect the dielec- 
tric properties shghtly (see the discussion below). Cell 
sizes (80 atoms for ab initio dynamics, 80 to 320 atoms 
for empirical) do not influence the structure appreciably. 

After cooling to K by damped dynamics, we fur- 
ther relaxed the structure using the ab initio Espresso^^ 
code (finding, incidentally, no major changes in structure 
or density). On the final structure, the Espresso code 
was used to perform linear-response density-functional 
perturbation theoryi^ calculations of electronic dielec- 
tric constant, effective charges and k=0 phonon modes, 
building up the dielectric constant (Eq. 1). All DFT cal- 
culations are done within the generalized-gradient ap- 
proximation (for details see Hflsl ). and will be indicated 
henceforth by the acronym GGA. 

III. RESULTS AND DISCUSSION 

A. Structural properties 

As no experimental structural data are available for 
the amorphous materials under consideration, a sanity 
check of our generated structures is to compare their 
properties with those of the crystal phase. The amor- 
phous samples appear to inherit many structural charac- 
teristics of the crystal i^ji. Their main feature is the ap- 
proximate conservation of an octahedral A-0 backbone, 
though with somewhat different octahedra distortion and 
degrees of cation disorder in LaAlOs and DyScOs. The 
structure-characterizing correlation function of the sim- 
ulated amorphous samples are collected in FigHJ Pre- 
liminarly we note that the intermixing of cation sublat- 
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FIG. 1: Pair correlation functions g(r) (left column) and an- 
gle distributions (right column) for amorphous LaAlOs and 
DyScOa. The g(r) functions are: (a), A-O; (c), R-O; (e), R- 
A. The angle distributions are: (b), 0-A-O; (d), O-R-O; (f) 
R-A-R. The functions for a-DyScOa are offset vertically for 
clarity. 



tices (i.e Ar antisites and such) is ruled out: no R^A 
cation exchange is observed in any of the simulations, 
and the hypothetical AlLaOa perovskite is metallic, 40% 
less dense and widely higher in energy than LaAlOa. 

Let us discuss FiglT] We refer to crystal or amorphous 
phases by the prefixes "c-" or "a-". In ideal perovskite 
the R (La, Dy) cation is 12-fold, and the A (Al, Sc) cation 
is 6-fold coordinated. The main feature of both a-phases 
is the A-0 g(r) function peaking neatly at typical A-0 
distances (Fig[T^) with a resulting A-0 coordination of 
6 as in the c-phases. The typical 0-A-O (e.g. 0-Al- 
O) angles are very close 90 and 180° in the crystals; as 
seen in Fig[T]3, in a-LaAlOs the angle distribution is a 
broad peak centered at 90°, whereas in a-DyScOs the 
same function is bimodal with peaks at about 80 and 
100°. This is due to the amplification of a Jahn- Teller- 
like asymmetry of the octahedron already incipient in 
the crystal)^. The a-phases can thus be characterized 
synthetically as disordered networks of apex-sharing oc- 
tahedra. 

The R-cations La and Dy are both approximately 6- 
fold coordinated in the respective a-phases. This means a 
severe reduction for La, whose coordination in c-LaAlOs 
is close to idealj- but almost no change for Dy which al- 
ready had** a low coordination of 6 in c-DyScOs due to 
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the large orthorhombic rotations. The Dy-0 g(r) is also 
sharper than the La-0 (FigH];), suggesting a lesser degree 
of disorder of the Dy sublattice. The loss of coordination 
of La is also signaled by the 0-La-O angle distribution 
(FigllJi), whose peaks at 60 and 120° are quite broad- 
ened but well-defined, whereas those at 90° and 180° are 
both washed out. This behavior is not unexpected as the 
ideal-perovskite angles are those internal to a cuboctahe- 
dron (60, 90, 120, and 180°): the 60° is the most robust 
as it corresponds to O pairs belonging to one single oc- 
tahedron, and correlates with the 120°, as does the 90° 
with the 180° one. 

Finally, the R-A functions show significant remnants of 
the correlation of the interpenetrating cubic sublattices 
of A and R cations in the c-phase: the g(r) main peak 
(FigdJj) is fairly well defined at the typical R-A distances, 
and the R-A-R angle function (FigHf ) still has a bimodal 
structure around the ideal-perovskite angles 70 and 110°. 



B. Polar izability 

With reference to Eq. 1, we note that the dipoles z 
generated by atomic vibrations determine, along with 
the vibrational frequencies discussed below, the static 
dielectric constant (Eq.l). One key ingredient of these 
dipoles is the ionic dynamical polarizability, i.e. as the 
time-integrated current flowing upon unit displacement 
of a specific ion.^ An unbiased measure thereof is pro- 
vided by the effective dynamical-charge tensor of atom i, 
Z*^^=r2/e {dPa/ dup), i.e. the derivative of the total po- 
larization vector with respect to the vector displacement 
of atom ii^ A drop in dynamical charges implies therefore 
a reduction of the ionic dielectric activity (Eq.l). 

We compare the present GGA dynamical charges, elec- 
tronic dielectric constants, and total static dielectric con- 
stants with their crystal counterparts^ii in Table HI To 
present the data synthetically, we average the a-phase 
charge tensors over all atoms of a given species, and de- 
compose them thereafter into s, p, and d components (i.e. 
spherical, traceless symmetric, and traceless antisymmet- 
ric). The s component, accounting for about 98% and 
95% of the tensor norm for cations and oxygen respec- 
tively, is reported in Table [J for the amorphous case. In 
the same Table we also report the s component of the di- 
electric tensors (these are only very mildly non-isotropic) . 
As we now discuss, the dynamical charges results of Ta- 
ble [T] are quite coherent with the structural properties 
described above, and with the general properties of the 
cations involved; charges on O largely adjust to obey dy- 
namical charge-neutrality sum rules. 

In a-LaAlOa, Al has the same non-anomalous dynam- 
ical charge as in the crystal. This is a token of its having 
no covalent bonds to support a dynamical charge flow 
upon displacement^ in either case. Al also has a rather 
crystal-like environment with a coordination of 6 and av- 
erage angle 90° (FigH^ and[T]3). For La, Z* is instead 
severely reduced: the disruption of ionic-covalent charge 



TABLE I: Average effective-charge, electronic, and static di- 
electric constant tensors in amorphous DyScOs and LaAlOs, 
compared to crystal values. See text for details. 
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crystal 
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24.0 


amorphous 


3.63 


3.83 


-2.43 


4.5 


22.1 


LaAlOs 
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crystal 


2.96 


4.37 


-2.44 


4.8 


26.4 


amorphous 


2.96 


3.88 


-2.38 


4.5 


27.1 



flow is clearly due to the major loss of coordination (12 
to 6, see FiglTJ;) also signaled by the 0-La-O angle dis- 
tribution in Fig[TJl. 

In a-DyScOa the effect is reversed. Dy has essentially 
the same Z* as in the crystal, in accordance with its 
having a similarly well-defined coordination and distance 
to O (Fig[TJ;). Note also the 0-Dy-O angle distribution, 
FigdJi, quite less blurred than the 0-La-O of a-LaAlOs. 
The A cation, Sc, suffers instead from an appreciable Z* 
reduction. As we discussed in Ref.|3, Sc has a sizably 
anomalous effective charge in the crystal, as a results of 
its employing d states in bonding to O. While the Sc co- 
ordination is still 6 in a-DyScOs, the covalent component 
of the dynamical polarizability is reduced by the strong 
octahedra distortion. We remind that the latter shows up 
in the bimodal angle distribution at 80 and 100° (Figllla) 
instead of a unimodal around 90°. 

In closing this section we recall our earlier reports^i^i^ 
of a sizable reduction of DFT effective-charge anomalies 
upon use of self-interaction correctionsJ^ This reduction 
lead to total dielectric constants smaller by ^10-20%, 
the agreement with experiment being typically improved. 
Due to the size of the cells involved here, we did not eval- 
uate the self-interaction-corrected charges in the present 
cases. The k's reported here may thus be expected to 
be somewhat larger than experiment. Of course, the 
amorphous-crystal comparison is unaffected. 



C. Dielectric and vibrational properties 

In TableUwe reported the average electronic and static 
dielectric constants. The a-phase static values are clearly 
in the vicinity of the crystal values, as we originally pur- 
ported to show. The electronic component is only slightly 
affected by disorder. Indeed, we find that both a-phases 
have a crystal-like gap in the electronic density of states, 
as expected in view of their strong ionicity and of the 
anion-like and cation-like character of, respectively, va- 
lence and conduction states. The ionic component, there- 
fore, is also approximately the same in both phases. To 
understand why this occurs despite the reduced effective 
charges, we analyze the frequency-dependent dielectric 
intensity, i.e. the individual (dimensionless) terms in the 
first equality of Eq.l. 

In Figl^lwe report the dielectric intensity (left panels) 



0.5 



0.5 





1 ' 1 ' 1 ' 1 ' 1 ' - 

(a) - 


■ (b) 


1 1 1 1 1 
IR cross-section " 








DyScO^ . 




Dielectric intensity — 
DyScO, 

III 

1 II L III 1 1 1 




IR 


- 
- 

" 1 


IR 

_ IVS/Ii / 

\ 


\ 




(c) J 


■ (d) 


IR cross-section ■ 


- 


Dielectric intensity 


- 


A LaAlOj . 




LaAlOj : 










/ IR 






^i^^ = 


R 








, 1 , 1 



100 200 300 400 500 
Frequency (cm ') 



200 



400 



600 



Frequency ( cm ) 



FIG. 2: Dielectric intensity (DI, left column) and IR spectra 
(right column) for a-LaAlOa and a-DyScOs. Panel a: DyScOa 
DI; b: DyScOa IR; c: LaAlOa DI; d: LaAlOa IR. In pan- 
els b and d the approximate positions of crystal vibrational 
modea^i^ is indicated. Labels stand for R=Raman, S=silent, 
Ii=infrared inefficient, IR=infrared. The dielectric intensity 
is dimensionless (the individual terms in the first equality of 
Eq.l). The IR spectrum, which has a 10-cm~^ lorentzian 
broadening built in, is in arbitrary units (same scale for pan- 
els b and d). 



and IR spectrum (right panels) of a-DyScOa (top pan- 
els) and a-LaAlOs (bottom panels). We indicated the 
approximate position of the calculated Raman, IR, silent 
modes of the crystals as horizontal bars superposed on 
the IR spectra of the a-phase. A detailed list is given in 
RefsH and 0- In Figl^l we report the dipole amplitudes 
decomposed by atomic species, i.e. the vectors z\ as de- 
fined in Eq.l, but obtained summing over one species 
only. 

Considering the IR spectrum which, we recall, is a mea- 
surable quantity, we note that the IR activity at low fre- 
quency is small. However, it is amplified in the dielectric 
intensity due to the l/w^ factor (Eq.l), so we concentrate 
on the low frequency region. To anticipate our conclu- 
sion, the dielectric enhancement at low frequency is due 
to IR-inactive (Raman, silent) or inefficient IR crystal 
modes that acquire IR character or enhance it due to 
disorder. (Inefficient IR modes in the crystal are those 
with weak intensity, say 3 orders of magnitude lower, 
than typical IR- active crystal modes - see Ref.|3.) 

Aside from minor differences of detail, both c-LaAlOs 
and c-DyScOa have their dominant IR modes in the range 
150-200 cm~^. As discussed previouslyj^ii these modes 
derive from the octahedron- vs-R sublattice vibration typ- 
ical of perovskites. Further IR modes around 300 to 
400 cm^^ and above - relatively uninteresting for us, 
being suppressed by the 1 /w^ factor in Eq. 1 - are intra- 
octahedron vibrations involving only O and the A cation. 
These features largely carry over to the a-phases, as con- 
firmed by the IR spectra (Fig[2h and d) showing peaks in 




200 400 600 



200 400 600 
Frequency (cm ) 

FIG. 3: (Color on-line) Average atom-specific dipoles in a- 
LaAlOa (panel a) and a-DyScOs (panel b), i.e. the effective 
mode charges zx (Eq.l) built up as a sum over each atomic 
species, normalized to the number of atoms in that species. 
Black curve: A cations; red (grey): R cation; green (light 
grey): O. 



approximate correspondence with crystal IR activity re- 
gions. The species-projected dipoles (Figl3]) also confirm 
this picture : in the 150-200 cm~^ region the dipoles 
generated, respectively, by the R cation and the A-0 
complex, are similar as expected. At higher frequency, 
instead, the displacements are exclusively A and/or O- 
like, and have intra-octahedron character. 

Let us discuss the low-frequency dielectric intensity- 
boosting modes starting with LaAlOs. We recall that in 
c-LaAlOa there exist three Raman modes (per primitive 
cell) related to the cubic-rhombohedral instability,'^ at 35 
(doubly degenerate) and 130 cm^^. They are character- 
ized by octahedra rotations accompanied by non-dipolar 
R-cation vibrations. In the disordered environment these 
modes mix with the IR counter-phase R-(A-O) motions, 
producing new vibrations with non-zero IR activity. The 
energies of these mixed modes start at about 40 cm"-'^, 
and boost the low-frequency dielectric intensity. The 
mixed modes can be seen as vibrations of the disordered 
R sublattice against A-0 units wobbling in a disorderly 
fashion about their center. This agrees with the dipole at 
low frequency being mostly contributed by the R cation, 
as shown in Figl3^. 

In c-DyScOs, the relevant modes^ are Raman at 100 
cm~^, silent at 70 cm~^, and inefficient IR also near 100 
cm~^. As we have discussed previously in Ref.Q - es- 
pecially in connection with Fig. 3 thereof - these three 
modes have no dipole (or just a tiny one) thanks to a 
delicate symmetry-induced compensation of different dis- 
placements. As such, they have a potential for develop- 
ing some dipole moment in a disordered environment dis- 
rupting the symmetric vibrational pattern. We therefore 
attribute the dielectric intensity boost at low frequency in 
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DyScOa to the activations of IR activity in these modes. 
Dy motions are again quite important in this energy re- 
gion. A peculiarity of DyScOa is that Sc motions (FiglSb) 
die out at frequencies above 300-400 cm"-'^ (unHke those 
of Al, Fig[3^). We note in passing that this is probably 
a characteristic feature of Sc-0 vibrations, as the dipoles 
are closely analogous to our findings for a-Sc203 to be 
presented elsewhere. 

We finally mention the results for the a-LaAlOa struc- 
ture generated by ab initio dynamics. The static dielec- 
tric constant for this structure is 22.4. While still in 
the vicinity of the crystal value, this result is apprecia- 
bly lower than for the empirical-potential sample (~27). 
The decrease is due to a (moderate) reduction in the 
ionic dielectric intensity at low frequency (the electronic 
part is unchanged, at 4.5). The reduction should be at- 
tributed to the lower average coordination of Al sites by 
about 10% (5.4 vs 6) in the ab initio sample compared to 
the empirical-potential sample discussed so. We suspect 
this to be an artifact caused by the much higher cool- 
ing rate of the ab initio sample (500 K/psec) compared 
to the empirical-potential one (50 K/psec) discussed in 
detail so far. (We recall, anyway, that all samples were 
optimized in structure at K after cooling.) 



Summary and acknowledgments 



In summary, we have shown that the static dielectric 
response in amorphous high-K RAO3 oxides is similar to 
that of their crystal phase, due to a combination of mod- 
erate polarizability reduction and of disorder-induced IR 
activation of non-polar low-energy modes, as well as to 
the inheritance by the a-phase of several c-phase struc- 
tural and vibrational features. We conclude that the 
dielectric-constant conservation upon amorphization of 
these materials is theoretically plausible. Given their 
very similar structural and vibrational properties at low 
frequencyi^, we suggest that our conclusions may quali- 
tatively apply to perovskites with cations from Group II 
and IV. 
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